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Abstract 

^' 

Cl^< We justify the validity of the discrete nonlinear Schrodinger equation for the tight-binding approx- 

imation in the context of the Gross-Pitaevskii equation with a periodic potential. Our construction 
of the periodic potential and the associated Wannier functions is based on the previous work [7], 
while our analysis involving energy estimates and Gronwall's inequality addresses time-dependent 
localized solutions on large but finite time intervals. 



/f ■ 1 Introduction 

'NT 

(N 



We consider the Gross-Pitaevskii (GP) equation with a periodic potential in the form 

i^t = -(t>xx + V{x)(t> + a\(t>\^4>, (1.1) 



where the solution (/> : M x M+ i-^ C decays to zero sufficiently fast as |x| — > oo, the potential y : R h-> M 
is a bounded 27r-periodic function, and the parameter a = ±1 is normalized for the cubic nonlinear 
term. In particular, we consider the piecewise-constant potential in the form 

H; t//^a _ / ^~^ X E (0, a) mod(27r) 

'^- yi^)-[ 0, xG(a,27r)mod(27r) ^^'^^ 

for some fixed < a < 27r and small e > 0. The asymptotic limit of small e represents the so-called 
tight-binding approximation, for which the potential V{x) is a periodic sequence of large walls of a non- 
zero width and the lowest bands in the spectrum of the linear operator L = —d'^ + V{x) are exponentially 
narrow with respect to e. According to the tight-binding approximation jl], time-dependent solutions 
of the GP equation (jl.ip are approximated by the time-dependent solutions of the discrete nonlinear 
Schrodinger (DNLS) equation in the form 

i<Pn = OL ((/"n+l + 0n-l) + (^ ii\<^nf '^n, (1-3) 

where a and /? are e-independent constants and the sequence {<i>n(S)\n&L represents a small-amplitude 
solution ^(x, t) evaluated at the periodic sequence of potential wells. 

We proved in the previous work [7J that stationary localized solutions of the GP equation in the form 
4>(x^ t) = $(x)e~*'^* with <1> G /7^(M) and lo ^ o"(L) are approximated for small values of e by stationary 



localized solutions of the DNLS equation in the form (j)n{t) = f^^e"*^*, where ^ € /^(Z) and Q is related 
to the rescaled parameter cv. Here and henceworth, we use the standard notations for the Sobolev space 
-ff^(M) of scalar complex- valued functions equipped with the squared norm 



{\cP\x)\' + \cP{x)f)dx 



and the space /^(Z) of vectors representing complex-valued sequences equipped with the norm ||0||/i(a) = 
Z^riez l'^"!- ^^ ^^^^ work, we extend our analysis to time-dependent localized solutions of these equations 
and prove that the formal tight-binding approximation of [1] can be justified for small values of e on 
large but finite time intervals. It will be clear from our analysis that the appropriate space for the time- 
dependent localized solution of the GP equation (jl.ip is associated with the quadratic form generated 
by operator —9^ -|- V{x) + 1. We denote this space by H^{R) and equip it with the squared norm 

UfniiR) = I {\<P'i^)\' + Vixm^)\' + [0(x)n dx. (1.4) 

Since V{x) > for all x € M, it is clear that ||(?!'|[hi(r) ^ II</'I|-hi(r)- 

Our analysis is closely related to the recent works on justifications of nonlinear evolution equations for 
pulses that exist in space-periodic potentials near edges of spectral bands [3] and in narrow band gaps 
of one-dimensional [TO] and two-dimensional [3] potentials. A similar work in the context of a nonlinear 
heat equation with a periodic diffusive term was developed in [9j with the invariant manifold reductions. 
Although the justification of lattice equations for the time-dependent solutions of dissipative (reaction- 
diffusion) equations can be extended globally for t > 0, the justification of the DNLS equation can only 
be carried out for finite time intervals because the GP equation is a conservative (Hamiltonian) system. 
Reductions to the DNLS equation on a finite lattice for a finite time interval were also discussed in the 
context of the GP equation with a A^-well trapping potential [2j . 

Methods of our analysis follow closely to arguments from [5j and rely on the Wannier function 
decomposition from [7\ as well as on energy estimates and Gronwall's inequality. The Wannier function 
decomposition is reviewed in Section 2. The energy estimates and the bounds on the remainder terms 
are studied in Section 3. The main theorem is formulated in Section 2 and proved in Section 3. 

2 Wannier function decomposition 

Let ui{x; k) be a Bloch function of the operator L = —d'^ + V{x) for the eigenvalue uJi{k), such that / € N, 
A; € T = [—2? 2] i^od(l), Ui{x + 27r; k) = ui{x; k)e^'^'"^ for all x € M, and the following orthogonality and 
normalization conditions are met 

ui,{x,k')ui{x,k)dx = 5i^v5{k-k'), V/,/' E N, VA;, A;' G T, (2.1) 

where di^ii is the Kronecker symbol and 5{k) is the Dirac delta function in the sense of distributions. To 
normalize uniquely the phase factors of the Bloch functions [6j, we assume that ui{x] —k) = ui{x; k) is 
chosen as a Bloch function for uji{—k) = idi{k) = 0Ji{k). 

Since the band function uJi{k) and the Bloch function n/(x; k) are periodic with respect to /c € T for 
any / E N, we represent them by the Fourier series 

u:i[k) = Y,^i^ne'^^^''\ ui[x;k) = Y,nUxy^^''\ V/GN, VfceT, (2.2) 

neZ neZ 



where the coefficients satisfy the constraints 

i^l,n = ^l-n = 1^1 -n, U^nix) = Ul^nix), Vn £ Z, V/ G N, Vx E M (2.3) 

and 

Ui,,{x)=ui^n-i{x-'^T^)=Uifi{x-2TTn), VnGZ, V/ GN, VxGM. (2.4) 

The real- valued functions ui^n{^) ai^e referred to as the Wannier functions. The following two proposi- 
tions from [7] summarize properties of the band and Wannier functions for the potential V{x) given by 
()1.2p in the limit of small e > 0. 

Proposition 1 Let V he given by 1^1. 2\) . For any fixed I G N, there exist eo > and e-independent 
constants (qjCOq, Cn > 0, such that, for any e G [0,eo); ^^^ hand functions of the operator L = —d'^ + V{x) 
satisfy the properties: 

(i) (band separation) min inf liOmik) — iOi o\ > Cq, (2.5) 

VmeN\{i} VfceT . I -- ' 

(ii) (band houndness) l'^/,ol ^ "^0) (2-6) 

{Hi) (tight-hinding approximation) p;_„| < c„e"e~~, n G N. (2-7) 

Proposition 2 Let V be given by hl.S^) . For any fixed I G N, there exists Cq > and e-independent 
constants Uo,Co,Cn > 0, such that, for any e G [0, eo); ^/^e Wannier functions of the operator L = 
—d'^ + V{x) satisfy the properties: 

(i) (houndness of norms) \\ui^n\\n^{R) I^Uq, n G N, (2-8) 

{ii) (compact support) \ui^q{x) — uq{x)\ < C^e, Vx G [0,27r], (2-9) 

{Hi) (exponential decay) \uifi{x)\ < CnS^e~~ , (2.10) 

Vx G [-27rn, -27r(n - 1)] U [27rn, 27r(n + 1)], n G N, 

where 

f 0, VxG [0,a], 

1 7fcsi^4;Fr^' VxG[a,27r]. 

Figure [1] illustrates the spectrum of L and the Wannier functions for the potential V in ()1.2p with 
a = TT and e = 0.5. The left panel shows the first spectral bands of L computed from the trace of 
the monodromy matrix [7j. The right panel shows the Wannier function ui^q{x) computed by using 
the integral representation iLifi{x) = Jj ui{x; k)dk and the finite-difference approximation of the Bloch 
function ui{x;k). The solid lines for the Wannier function 'Uio(x) approach the dotted line for the 
asymptotic approximation (j2.1ip with / = 1 as e gets smaller. 

We use the set of Wannier functions {iii^n}n£Z for a fixed / G N to represent formally solutions of the 
GP equation (jl.ip in the form 

cb{x,t)=n^/^{ipo{x,T)+nip{x,t))e-'^''°\ T = fj,t, /i = ee-t, (2.12) 

where 

ipo{x, T) = Y, <Pn{T)ui^n{x). (2.13) 

nez 
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Figure 1: Left: the band-gap structure of the spectrum (y{V) for e = 0.5. Right: the Wannier functions 
u\id{x) for e = 0.5 and e = 0.25. The dotted Hne shows the asymptotic approximation (|2.1ip . 



Let us assume that the sequence {(/'ninez satisfies the DNLS equation ()1.3p with a = ^^ and /3 

|4 



|ii;^o|l24(-iR)- Since the Wannier functions satisfy the ODE system from [7] 



(2.14) 



n'e2 



we obtain an inhomogeneous PDE system for the function {p{x,t) in the form 

lift = -if XX + V{x)if - LOiflif + - X^ X^ <^l,m {(t>n+m + 4>n-m) Ul,n 



fJ' 



neZm>2 



+ 0- I \<po + Hifl'^ {^0 +^(/7) -/3^ |0n| 
V neZ 



'0nWL- 



(2.15) 



The term Ic/jqPv'o gives projections both to the selected l-th. spectral band and to its complement in 
L'^{M.). The following two lemmas allow us to control both projections. 

Lemma 1 Let Ei be the invariant closed subspace o/L^(M) associated with the l-th spectral band and 
assume that Ei n E^ = for a fixed I € N and all m ^ I. Then, {un,i,Un',i) = ^n,n' for any n,n' £ Z 
and there exists constants rji > and Ci > 0, such that 



ui 



,n(x)| < Qe-'''l^-2'^"l, Vn G Z, Vx G M. 



(2.16) 



Moreover, if cf) £ l^C^), ui,n ^ TC^(M.), and (j){x) = Xlnez '?^n^^n(^) f'^''" ^ fixed I G N, then cp G Ei, 
{(t),i;)=0,yijeUm^iEm,and^en\R). 

Proof. The orthogonality and exponential decay of Wannier functions follows from the orthogonality 
relations (j2.ip and complex integration (see [7] for the proof). The assertion that (p € Ei and {(j), ip) = 0, 
V^ G Um=/=iEm follows from the L^ spectral theory for the operator L = —d"^ + V{x) (if cf) G /^(Z), then 
(f) G /■^(Z) and (j) G L^(M)). The assertion that (j) S 'hO'^M) follows from the triangular inequality. D 



Remark 1 According to property (i) of Proposition [U the l-th spectral band for a fixed I G N becomes 
disjoint from the rest of the spectrum of L in the asymptotic limit e — > 0. According to property (i) of 
Proposition [21 ui^n ^ 'H^(M) for all n G N uniformly in e > 0. Therefore, the assumptions of Lemma [T] 
are satisfied for sufficiently small e > 0. 

Lemma 2 Let 11 be an orthogonal projection from L^(IR) to Ei C L^(IR). There exists a unique solution 
if ^Ti} (M) of the inhomogeneous equation 

i-dl + Vix) - ui,o) ^ = {I-U)f, (2.17) 

for any f G L^(M), uniformly in s >0, such that (v?, ^) = 0, ^tp G Ei. 

Proof. By property (i) of Proposition [H if / G L^(M), then ip G L^(M) uniformly in e > 0. By property 
(ii) of Proposition [H we obtain 

WWIhr) + ll^'/Vlli2(M) < |(2'/,olll'^lli2(M) + \ivJ)\ < cil/lli2(M), (2.18) 

where the constant C > is e- independent. Therefore, if / G L^(M), then ip G 7Y^(M) uniformly in 
e > 0. Uniqueness of 99 follows from the fact that the operator L — uJifi is invertible in L^(]R)\S;. D 

We can also use the following elementary result. 
Lemma 3 The space 7^^(M) forms Banach algebra under the pointwise multiplication, such that 

\/u,v G n^im) : ||™||7^i(K) < C\\u\\'Hi(R)\\v\\'n^R), (2.19) 

for some C > 0. 

Proof. The result follows from the representation HuH^imN = II^II|/i(k) + 11^^ ^^IliafR) ^iid the Sobolev 
embedding theorem ||u||^oo(ig) < C||n||j|^i(]g') for some C > 0. D 

Let us return back to the evolution problem (j2.15p and decompose the solution f{x, t) into two 
parts (p{x,t) = ipi{x,T) + ip^Xjt), where ipi is a solution of the inhomogeneous equation (j2.17p with 
/ = —a\ipQ{x,T)\'^ipQ{x,T), while ip satisfies the evolution problem in the abstract form 



where 



#i = (L - c:;/,o) i^ + /^^(0) + f^(TN{cf>, -0), (2.20) 



R{$) = —^^ '^l,m i4>n+m + (t>n-m) Ul,n + - U\(fo\'^(fo " /? ^ \(pnf(pnUl,n (2.21) 

^ neZm>2 ^ \ neZ / 



and 



+fi {2\ipi + VI Vo + (<^i + i'f^o) + ^^If/^i + ^\\vi + ip), (2.22) 

with (fo = J2n&N 'i^nUi,n, Vi = -(^{^ -'^)iL - iVi^y^ (I - U)\ipQ\'^ipQ, and a = ±1. The following lemma 
gives a bound on the vector field of the evolution problem (j2.20p . 



Lemma 4 Let Ds^ C /^(Z) be a ball of finite radius 5i centered at G l^{'^), Ds^ C H^(R) be a ball of 
finite radius 62 centered at G Ti.^(W) and -R^o cM. be an interval of small radius fj,o centered ai G M. 
Then, for any /i € (0, /Uq), ||</>|I«i(z) £ [0,5i) and ||V'llwi(R) ^ [0)'^2)) there exists ^-independent constants 
Cr, Cn > such that 

(2.23) 



ll^(0)llwi(R) < CR|!0|bi(z), ||7V(0,V')|lwi(R) < C^ ||0lbi(Z) + ll^llwi 



Proof. By the last assertion of LemmalU if ^ € l^{7,) and (p{x) = J2n€Z 4'n'i''i,n{x) for a fixed / G N, then 
(j) G 7Y^(M). Therefore, there exists Cq > such that ||'/3o||7^i(ir) < Co||0||/i(2). By Lemmas [2] and [3l there 
exists Ci > such that ||¥'i|l7^i(]R) < Ci||0||^i(.^w where we have used the fact that |||</'op¥'o||wi(M) ^ 
C^||(^o||f^i(K) for some C > 0. The vector field -R(<^) can be represented by -R(0) = J2n&z''^n{<P)ui,n{x) , 



where 






m>2 



(7 

'^Z,m Wn+m + (pn-m) H 

A' 



E 



i^, 



n,ni ,712 iJT-a r"ni yn2 Y'na ; 



(ni,n2,n3)eZ3\{(n,n,n)} 



where Kn,nx,n2,n3 = {ui,n,ui,niUi,n2^i,n3)- The first bound in (j2.23p is proved if r E 1^{Z) for every 
cf) E /""^(Z) and the map r{cj)) is uniformly bounded for small /i > 0. The first term in r(0) is estimated 
as follows 



/ ^ i^/^m (<;/'n+m + '/'n- 



m>2 



< 






/ ^ / ^ \'^l,m+n\\(t>n\ < Ki\\<t)\\i 

^l.^^ neZmeZ\{0,l,-l} 



where Ki = sup„g2\ ro^i^_i} X^mez 1'^^,"+™!- Since u)i{k) is analytically extended along the Riemann 
surface on /c E T (by Theorem Xin.95 on p. 301 in [8j), we have loi E H^{T) for any s > 0, such that 
Ki < 00. The second term in r(0) is estimated as follows 



/ ^ ^n,ni,n2,n3,Yni^'n2VTi2, 

(ni,n2,n3)eZ3\{(n,n,n)} 



< 



\Kn 



7 ^ 7 ^ i-"-n,ni,n2,n3lirni ||S^n2ll¥^n3l 

^1(2;) neZ (rai,n2,n3)eZ3\{(n,n,n)} 



< i^9 



where K2 = sup(„^^„2^„3)g2\|(„ „_„)} X^^^^ |-?^n,ni,n2,n3l- Using the exponential decay (|2.16p, we obtain 

j;|n,,„(x)|<Cz^e-'''l^-2-"l<A 

neZ nGZ 

for some Ai > {) uniformly in x E M and 

^ |i^n,ni,n2,n3l < A / |^«,ni (2;) | |u/,n2 (2^) I |^tZ,n3 (^) l^^; < A||lt/,ol 



L4(iR,)||ni,o||L2( 



uniformly in {ni,n2,n^) E Z^. By the Sobolev embedding theorem and property (i) of Proposition [21 
ll^i,ollL4(R) ^ C'||{t; o||_ffi(R) < C'll^z,oll'Hi(iK) fo'^ some C > 0, such that K2 < 00. Therefore, the norm 
11-^(0) li'Hi(R) is bounded from above by the norm ||0||ii(z)- 

To show that the constant Cr is uniform for small /i > 0, we use Propositions [T] and [2J By property 
(iii) of Proposition [H u^i^m = 0(^'") for all m > 2, such that Ki/ jj? is uniformly bounded for small ^. 



By property (iii) of Proposition [21 K„,ni,n2,n3 = O (;/l"i-«l+l'^2-"l+l"3-n|+|n2-ni|+|n3-ni|+|n3-n2|^ for ^n 
ni,n2,n3 E Z^\{(n, n,n)}, such that K2/H, is uniformly bounded for smah y^. Thus, the first bound in 
<n:7m is proved. 

The second bound in ()2.23p follows from the fact that both 7^^(M) and /^(Z) form Banach algebras 
with respect to pointwise multiplication. As a result, if G ^^(^) and <^(T) is a solution of the DNLS 
equation ([TS]), then dx^) G /H^) and if ipo,ipi G W^(R) and ^ G Z^(Z), then A^(0,'0) maps V G H^(]R) 
to an element of rO-{W). D 

We can now prove that the initial-value problem for the time-evolution equation (j2.20p and the 
initial-value problem for the DNLS equation (|1.3p are locally well-posed. 

Theorem 1 Let (t>{T) G C^(1R, /^(Z)) and ipQ G 'H^(M). Then, there exists a io > and a unique 
solution 'ip{t) G C^{[0,tQ],'H^{M.)) of the time- evolution problem 112. 20\) with initial data ip{x,0) = iPq{x). 

Proof. Since L is a self-adjoint operator, the operator e~^^^^~'^'-'^' forms a strongly continuous semi-group 

and 

|[g-it(L~i,,o)|| <^^^ 

for some Kq > uniformly in t G M+. Using the variation of constant formula, we rewrite the time- 
evolution problem ()2.20p in the integral form 

^(t) = e-**(^-'^''0)Vo + I e-*(*-")(-^-'^'>o) LR{$is)) + ^iCjiV(0(s), -0(5))) ds. (2.24) 

By using bounds (j2.23p on R{(f)) and N{(f), ip) and the contraction mapping principle for sufficiently small 
to > 0, one can show with a standard analysis that there exists a unique solution ^(t) G C"'^([0, tQ\,li}(^)) 
of the integral equation (j2.24p . D 

Theorem 2 Let (pQ G 1^{Z,). Then, there exist a Tq > and a unique solution (t){T) G C^([0, Tq], /^(Z)) 
of the DNLS equation ( tJ.3)) with initial data 0(0) = cj>q. 

Proof. By the variation of constant formula, we have 



${T) = $^-i f (^aA${s) + apT{${s))) ds, 



where (A0)„ = (pn+i + (pn-i and (r(0))„ = \4>n\'^4'n- Since 1^{1j) forms a Banach algebra, the right- 
hand-side of the integral equation maps an element of /^(Z) to an element of /^(Z). Therefore, there 
exists a unique solution (p{T) G C^([0, Tq], /^(Z)) of the integral equation for sufficiently small Tq > 0. 
D 

We can now formulate the main theorem of our article. 

Theorem 3 Fix i G N and let ${T) G C^{[0,To],l^{Z)) be a solution of the DNLS equation < f73|) with 
initial data 0(0) = 4>q satisfying the bound 



4>o -/i^^^^0n(O)-u/,„(a 



< Con^^^ (2.25) 

-Hi(M) 



for some Cq > 0. Then, for any fi G (0, ^Uq) with sufficiently small fj,Q > 0, there exists a ^-independent 
constant C > such that equation ( [J.i|) has a solution (f){t) € C^([0, To//i],'H-'^(M)) satisfying the hound 



VtG[0,ro/;u]: 






< C7/i3/2. 



(2.26) 



Wi 



Remark 2 Since // = ee "^'^j the finite interval [0, 7b//x] is exponentially large with respect to parameter 
e, similar to the bounds obtained in [2j. 

Theorem [3] is proved in the following section. 



3 Bounds on the remainder terms 



We develop the proof of Theorem [3] by using the energy estimates for the time-evolution problem (j2.20p 
and Gronwall's inequality for a scalar first-order differential equation. The GP equation (jl.ip has two 
conserved quantities 



Q(<^) 



^dx, E{(t)) 



j^[\<P,\^ + V{xm^ + \a\ct>\''^dx, 



(3.1) 



which have the meaning of the charge and energy invariants, such that Q{(t)) = Q{(po) ^-^d E{(j)) = E{(J)q) 
for any solution (j){x,t) starting from the initial data (j)o{x). We shall consider the quantity Eq{'iP) = 
llV'll'Hi(R)) which is not a constant in t if ip{x,t) satisfies the time-evolution problem (j2.20p . The time 
evolution of Eq^ip) obeys the following estimate. 

Lemma 5 Let 4>{T) € C^(R4-, /^(Z)) he any function and ip{t) G C^{[0,to],TC^{]R.)) he a local solution 
of the time- evolution prohlem h2. 20\) for some to > 0. Then, for any /i E [0, 1] and every M > 0, there 
exist a ^-independent constant Ce > such that 



|^.(^) 



<liCEm\\iHz)+EQ{i^) 



(3.2) 



as long as Eq{'iP) < M . 



Proof. By direct differentiation, for any i/^(t) G C^([0, to])'^^(J^))) we have 
d_ 



in / il)xRx{(l}) - i)xRx{(t>) dx - i^a 
ifM f {1 + V{x)) {i)R{$) - ^R{$)\ dx 

1 + V{x)) Un{$, ^) - ipN{$, iP)] dx 



ipxNxicf), ip) - ipxNxi^, ip) ) dx 



-ifia 



Using the Cauchy-Schwartz inequalities and the bounds (j2.23p of Lemma HI we obtain 
d 



dV 



\2 



< 6^11^1 



H^ 



{CR + CN)U\\n{z)+CN\ 



llwi 



where o" = ±1 has been used. By canceling one power of |!V'II'K1(r)> we arrive to the bound (j3.2p . D 
Theorem [3] is then a direct consequence of the following corollary. 



Corollary 1 A local solution ip{t) G C^{[0,TQ/fj],T-i^{]R.)) of the time- evolution problem 112. 2U\) for any 
V'(O) G ^^(M) and any ${T) G CH[0,To], ^^Z)) satisfies the hound 

sup ||V'(t)||7^i(R)< [||V'(0)||wi(R)+CEro sup ||0(T)||,i(2))e^^^«. (3.3) 



Proof. By using the bound (13. 2p . we obtain 

i^gCV'W) < i^Q(V'o) + A^Cb y (||0(/"S) 11^1(2;) + EQ{ij{s))) ds. 
The bound (j3.3p follows by Gronwall's inequality on t G [0,To//i]. D 
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